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Abstract. The escape rate of cosmic ray (CR) particles from radio galaxy cocoons is a problem of high astrophysical relevance: 
e.g. if CR electrons are stored for long times in the dilute relativistic medium filling the radio cocoons (radio plasma in the 
following) they are protected against Coulomb losses and thus are able to produce a significant non-thermal Comptonisation 
signature on the CMB. On the other hand, CR protons and positrons which leak out of radio plasma can interact with the 
ambient medium, leading to characteristic gamma ray radiation by pion decay and pair annihilation. In order better understand 
such problems a model for the escape of CR particles from radio galaxy cocoons is presented here. It is assumed that the radio 
cocoon is poorly magnetically connected to the environment. An extreme case of this kind is an insulating boundary layer of 
magnetic fields, which can efficiently suppress particle escape. More likely, magnetic field lines are less organised and allow the 
transport of CR particles from the source interior to the surface region. For such a scenario two transport regimes are analysed: 
diffusion of particles along inter-phase magnetic flux tubes (leaving the cocoon) and cross field transport of particles in flux 
tubes touching the cocoon surface. The cross field diffusion is likely the dominate escape path, unless a significant fraction of 
the surface is magnetically connected to the environment. Major cluster merger should strongly enhance the particle escape by 
two complementary mechanisms, i) The merger shock waves shred radio cocoons into filamentary structures, allowing the CRs 
to easily reach the radio plasma boundary due to the changed morphology, ii) Also efficient particle losses can be expected for 
radio cocoons not compressed in shock waves. There, for a short period after the sudden injection of large scale turbulence, the 
(anomalous) cross field diffusion can be enhanced by several orders of magnitude. This lasts until the turbulent energy cascade 
has reached the microscopic scales, which determine the value of the microscopic diffusion coefficients. 
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1. Introduction 

1.1. Motivation 

The outflows of radio galaxies fill large regions in the inter¬ 
galactic space with relativistic, magnetised plasma. The fate 
of this radio plasma is unclear, since it rapidly becomes unde¬ 
tectable for radio telescopes due to the radiative energy losses 
of the higher energy electrons, which emit the observable syn¬ 
chrotron emission at radio bands. Lower energy electrons, and 
any possible present relativistic proton population, may reside 
for cosmological times, unless they are able to escape by spatial 
diffusion. 

The escape of CRs out of radio plasma is of high astro- 
physical relevance. If CR electrons are able to leak out of the 
cocoon, they are possible seed particles for the Mpc-sized ra¬ 
dio halos in clusters of galaxies if they are re-accelerated in 
the cluster turbulence to radio-observable energies (Giovannini 


line, if they would be able to leave the radio cocoon and inter¬ 
act with a dense ambient intra-cluster medium (Furlanetto and 


Loeb, 2002). 


On the other hand, if the relativistic electrons are efficiently 
confined for cosmological time-scales in radio cocoons, they 
are shielded from severe energy losses by Coulomb interac¬ 
tion with the environmental gas. In such a case they would 
be long-lived and would therefore be able to produce a non- 
thermal Co mptonisation signatu r e in the cosmic microwave 
backg round ( EnBlin and Kaiser , 2000 ; fcnfilin and Sunyaev , 
2002 ). Further, if old, invisible radio plasma with confined rela¬ 
tivistic electrons is dragged into a shock wave of the large-scale 
structure formation, its radio emission is possibly revived and 
forms the observed cluster radio relics ( EnBlin and Briiggen , 
2002, and references therein). 


The recent detections of several ghost cavities in galaxy 


et al., 1993; Brunetti et al., 2001). If radio cocoons release CR clusters (for lists of recent detections e.g. EnBlin and Heinz 


protons into the intra-cluster medium, secondary CR electrons 
are produced in hadronic interactions with the background gas, 
which also could explain the existence of cluster radio halos 


2002; Soker et al., 2002), which are often but not always radio- 


(Dennison, 1980). If radio plasma contains a significant frac¬ 


tion of positrons, they could lead to a detectable annihilation 


emitting, support the assumption that the environmental gas 
and the bulk of the relativistic plasma stays separated on a 
timescale of 100 Myr. However, leakage of higher energy par¬ 
ticles is not excluded by these observations. 
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CR escape from radio plasma can also play a role in our 
galaxy. Microquasars are proposed to contribute characteristic 
spectral features to the galactic CR spectrum, if CRs are able 
to leave the ejected plasma (Heinz and Sunyaev, 2002). 


1.2. A simplified approach 

The usual description of CR escape from systems like our 
Galaxy is done in terms of leaky box models, where the CR 
particles trapped in some region are attributed a characteristic 
escape time or frequency (t = v^. ape ). This escape time is of¬ 
ten empirically determined (e.g. for th e CR esc a pe fro m the 
Galaxy by radioactive CR clocks (e.g. Cesarsky, 1980)). It is 


gyro-radii scales. In any numerical example a Kolmogorov- 
like turbulence spectrum is adopted. The detailed depen¬ 
dency on the turbulence wave types (e.g. Michalek and 
Ostrowski 


1997, 1998; Michalek et al. 


1999; Chandran, 


the goal of this work to provide an estimate of this escape time 
and its dependency on the geometrical and plasma parameters 
of the system under investigation. The focus of this work is 
cocoons of radio galaxies, but the presented model may have 
other applications. 

In many applications, the physical parameters will not be 
sufficiently constrained to allow an accurate calculation of the 
escape time. Nevertheless, insight into the dependencies on pa¬ 
rameters such as turbulence, fraction of the magnetically open 
surface etc. allows statements of the relative CR escape rates in 
different situations to be made. This should help to formulate 
hypotheses about the conditions in which CR escape is effi¬ 
cient, which may be tested observationally. Examples of such 
tests were given in Sec. 0 

The following simplifications are used to compile the CR 
escape model: 


2000b) is put into a phenomenological “fudge factor’’. 

- The macroscopic cross field diffusion coefficient (/e a ) is as¬ 
sumed to exceed the microscopic one by orders of magni¬ 
tude due to magnetic field line wandering. This allows CRs 
that are rapidly diffusing along the field lines to be trans¬ 
ported across them (Rechester and Rosenbluth, 1978). This 
also ensures that the length scales of the diffusion problem 
are always shorter than length scales on which a magnetic 
flux tube loses its identity due to field line wandering. This 
should justify the use of the flux tube picture. 

1.3. The structure of the paper 

The structure of this article is the following: In Sect. [|the the¬ 
oretical tools for the description of the CR propagation in in¬ 
homogeneous media are compiled. The inter-phase CR trans¬ 
port is analysed in Sect. |j. Three escape routes are considered: 
the penetration of an isolating boundary layer (Sect. 3T), the 
parallel diffusion along inter-phase magnetic flux tubes (Sect. 
[T2| ), and the cross field escape (Sect. [k3|). In Sect. [| the CR 
escape from radio plasma cocoons is investigated and in Sect. 
U the main findings of the paper are listed. In Appendix A the 
spatially homogenised transport equation of CRs in a small- 
scale inhomogeneous medium is derived. Appendix B contains 
a glossary of the frequently used symbols and abbreviations. 


Convective CR transport, e.g. by plasmoids which detach 
by reconnection events from the source region, is neglected. 
The detailed magnetic structure of the source interior is not 
modelled and is treated as being homogeneous. Only the 
topological properties of the magnetic flux tubes touching 
the source surface are modelled. 

The magnetic field at the source region surface is virtually 
split into flux tubes with diameters of the perpendicular au¬ 
tocorrelation length of the magnetic field. The transport of 
particles along and across the flux tubes is described as 
a diffusion process. A similar phenomenological descrip- 
tion of CR transpo rt within flux tubes can be found in e.g. 
jChandran ( 2000a ), where limits of such a picture are also 
discussed. 

It is assumed that such a surface flux tube bends at some 
point into the interior of the source region. This allow CRs 
from the interior to enter the tube there and to follow it to 
the surface. In the case that the magnetic fields are more 
onion-like in structure, the particle transport will be much 
more suppressed than predicted with this model. A brief 
discussion of such a hypothetical situation is given in Sect. 


3.1 


- A simplified description of the diffusion process is 
adopted: The microscopic diffusion coefficients ( a :||, k ± ) are 
parametrised. The spectrum of magnetic field fluctuations, 
which strongly determines the transport coefficients, is as¬ 
sumed to be a single power law, connecting the length 
scales on which field line wandering happens, down to CR 


2. Cosmic ray diffusion 

2.1. Diffusion in inhomogeneous magnetic fields 

Since the magnetic field topology dominates the mobility of 
CR particles, adapted coordinates are chosen. The length x is 
measured along the (local) mean direction B of a field bundle 
with local field strength M = B + 6B. Although individual field 
lines may leave the flux tubes, due to magnetic fluctuations SB, 
this mean field flux tube is well defined and gives an appropriate 
local coordinate system. 

The analysis is restricted to particles with gyro-radii much 
smaller than the length-scale of the magnetic fields. In this 
case the cross field diffusion is orders of magnitude slower 
than the diffusion along the field lines. Particles can be re¬ 
garded as being confined in a flux tube. Since the diameter 
and the field strength of the tube can change as a function 
of x, it is convenient to work with the number of particles 
F per flux tube length dx and magnetic flux dtps = Bdydz, 
instead of the usual volume normalisation dV = dxdydz = 
dxdtpB/B. The particle phase space distribution function is 
F(x,p,p,t) = dN/(dx d(f>B dp dp), and its source density is 
Q(x,p,p,t ) = dN/(dx d(pB dp dp), p is the momentum of the 
particle, v its velocity, m is its mass, and c is the speed of light. 
p - p x /p is the cosine of the pitch angle between B (or x ) 
and p. The particle distribution is assumed to be rotationally 
symmetric with respect to B, therefore the azimuthal angle has 
been integrated out. The particles entering or leaving the flux 
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tube are included in the source term Q and in the loss term 
-F/t in the Fokker-Planck equation for F: 


dF d 

-* + aC" F) ' 


dp 


(P F) ■ 


d_ 

dp 


(A iF) 


d „ d „ F „ 
-p^Tp - +Q - 


trum of plasma waves (e.g. 

Fokipii, 

1966 

1967 

; Kulsrud 

and Pearce, 

1969;F 

[asselmann and Wibberenz 

197C 

; Skilling, 

1975b; Schlickeiser 

, 2002 

, and references in the latter). Further 




V) 


dB 
B dx 


Kit: pitch and Pluskir? (1995) and Chandran (2000a) investigate 


f(x,p,t ) 
q(x, p, t ) 


dN 


dxd(f>sdp 

dN 

dxdipsdp 


-jG 

-jG 


dp F{x,p,p, t ), 


dp Q(x,p,p, t). 


q 


This equation describes the diffusive transport of particles 
along the flux tube. In the case that small-scale variations of the 
coefficients in this equations exist, effective large-scale coeffi¬ 
cients can be estimated, as shown in Appendix A. The diffusion 
coefficient is 


( 1 ) 


K\\(X,P) 


v (p) 


is the Fokker-Planck pitch angle diffusion coefficient, 
which can in principle be calculated from quasi-linear the¬ 
ory of plasma wave-particle interaction for a given spec- 


ft 


2\2 


(i -p~y 


8 J-i D m (x,p,p) 
and the pitch angle averaged momentum losses are 


p(x, p) 


1 r 1 

2 Jdpp(x,p, 


p). 


( 6 ) 


(7) 


Fokker-Planck diffusion coefficients, which describe momen¬ 
tum changes of the particles, can be found in in the literature 
(e.g. [Schlickeiser , 1989a bj). Here, momentum diffusion is ig¬ 
nored due its small impact on spatial transport processes, and 
only continuous energy or momentum losses are included by 

P- 

The continuous pitch angle changes can be calculated from 
the adiabatic invariants of a particle with charge Z moving 
in spatially slowly varying, and temporally constant magnetic 
fields. Specifically, these are the linear momentum p and the 
magnetic flux within a gyro-radius <pB(r g ) = n r 2 ^ B, where 
r g (x,p,p) = (1 - p 2 ) l/2 pc/(ZeB(x)) is the gyro-radius. One 
gets 


In order to demonstrate that Eq. is the proper transport equa¬ 
tion the volume density of particles with momentum p is in¬ 
troduced: g(x,p,t) = dN/(dVdp ) = B(x) f(x,p,t). Eq. |]then 
transforms into 


9g , 9 ,. , R d (k\\ 3g\ g 
^ + ^ P8) = B MBa-x)~r + Bq - 


( 8 ) 


( 2 ) 


the weak scattering regime, in which the mean free path of a 
CR is large compared to typical length-scales of the inhomo¬ 
geneous magnetic fields and CR density scale-length. Here, a 
more conventional standpoint is adopted, by assuming that the 
particle distribution is rapidly isotropised by plasma wave in¬ 
teractions. |] Therefore the anisotropic part of the distribution 
function is small compared to the isotropic one. The pitch an¬ 
gle integrated distribution and injection densities are defined 
by 


This e quation is consistent with Eq. 2 of |Cesarsky and Volk 
( 1978 ). Diffusion of g vanishes for dg/dx = 0. This means that 
diffusion tries to approach a state where the space density of 
particles with identical momentum is constant, as it should do. 
Due to variations in the diameter of the flux tube this can result 
in a non-constant / as a function of position. 

2.2. Diffusion coefficients 


The parallel diffusion coefficient depends on the pitch an¬ 
gle diffusion coefficient according to Eq. |(| By defining the 
particle-wave scattering frequency v t . 


,,-i 


f 


1 

djj ■ 


2\2 


(1-/0 


D m (x,p,p) 

one can write the parallel diffusion coefficient as 

KBohmCr p) 


K\\<,X, p) 


s(x, p) 


(9) 


( 10 ) 


(3) 

(4) 


With the help of a quasi-linear approximation the evolution 
equation for / can be derived from Eq. [j] following the cal¬ 
culation steps described in Schlickeiser (1989a): 


Here, e(xp) is the ratio of scattering frequency v ; , to the 
gyro-frequency Q(x, p) - v(p)/r g (x, p,0), and KBohm(x,p) = 
v(p) r g (x, p, 0)/3 = v(p) p c/(3 Ze B(x )) is the Bohm diffusion 
coefficient. v fl can be regarded as the decay time of the paral¬ 
lel particle velocity autocorrelations, which in the following is 
assumed to be identical to the decay time of the perpendicular 
velocity autocorrelations. This latter assumption allows us to 
write the microscopic perpendicular diffusion coefficient as 


(5) 


<</ <v , . I w>\ / 

* + w> ip f) ‘ 

1 This is supported by the observation that the pitch angle diffu- 
sion coefficient is the fastest of all the Fokker-Planck coefficients (e.g. 


1 + E‘ 


' ^Bohm 


(ID 


(Bieber and Matthaeus, 1997). 


Schlickeiser, 2002). The relatively fast pitch angle scattering should 


help to maintain a nearly isotropic CR distribution - as observed in our 
own galaxy. A criterion to test if pitch angle scattering is indeed suf¬ 
ficient to establish an isotropic pitch-angle distribution - as assumed 
throughout this paper - is given later. 


The scattering frequency v^ depends on details of the un¬ 
derlying plasma turbulence on scales comparable to the gyro- 
radius of the particle. There are two contributions important 
to this small-scale turbulence: The first is the Kolmogorov (or 
Rraichnan) cascade of large-scale turbulent kinetic energy to 
smaller length-scale, and the second is turbulence induced by 


CR streaming (Wentzel, 1968, 1969; Skilling, 1975a). Since 
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the focus of this work is on the transport of poorly connected 
regions, the amount of CR streaming is expected to be low. In 
the following numerical examples, only external Kolmogorov- 
like turbulence is assumed. In other cases, the theory of this ar¬ 
ticle can still be applied if the appropriate diffusion coefficients 
are used. 

A simplified parameterisation of the scattering frequency is 
adopted here. It is assumed that 


parametrisation we get 


tip) 


so SB 2 (r g (p))/B 2 , 


SBr{l) = (J 0 B~ (l/I B ) y , 


sumed that a single power-law describes the magnetic fluctua¬ 
tion spectrum from the scales important for field line wandering 
down to the CR gyro-radii. 

The parameter 60 and eq play an important role for many 
of the addressed questions and their expected values should 
briefly be discussed here, although our knowledge of these 
quantities is still very limited. So is roughly speaking the ra¬ 
tio of the magnetic power on the largest length-scale, on which 
the power-law spectrum of the inertia range of the turbulence 
is valid, to the total magnetic energy density. For a sharply 
peaked, single power-law magnetic spectrum it is <5o ~ 1. 
For a broken power-law spectrum, or a very broad maximum 
above the inertia-range length-scales, we expect ho <s 1. eo 
gives the efficiency of CR scattering per magnetic power on 
length-scales of the order of the particle gyro-radius, thus the 
power on scales which can resonate with the particle gyro- 
orbit. For an Alfven wave spectrum with slab-like geometry 
(waves-vectors are mainly parallel to the magnetic main di¬ 
rection) one expects eq ~ 1 .[] However, the nature of MHD 
turbulence might be anisotropic on small scales in the sense 
that mainly wave-modes with wave-vectors perpendicular to 
the main field are populated ( Sridhar and Goldreich , 1994 ; 
Goldreich and Sridhart 1997 ). In that case a strong reduction 
of the scattering frequency ca n be expec t ed (gp <g 1) (a quasi- 
linear estimate is provided by Chandran, 2000b, however it is 


4catt(P) — 


h 

£o do 
10~ 3 kpc 
so do 


i-r 


(14) 


( pc F / Ib V I ZB 
\GeV/ \kpcj \pG 


( 12 ) 


This is sufficiently small in most of our cases to allow the dif¬ 
fusive approximation to be used. 

The resulting diffusion coefficients for CR particles are then 


where so is a fudge factor, which allows one to trace and correct 
the error made in this simplification, and 


*11 = 


pc 


vIb 


3 so do \Ze B Ib 


i-r 


(15) 


(13) 


= 3.2- 10 


.28 


cm 


is the power in the magnetic fluctuations on scale / (and 
smaller). In Kolmogorov-like turbulence y = 2/3, which will 
be used in numerical examples. The reference scale Ib is 
the largest length scale which contains significant magnetic 
power (compared to the power-law given by Eq. |l3|), so that 
6B 2 (Ib) = {SB 2 ) is the total magnetic fluctuation power. This 
length scale is also of the order of the coherence length of 
the fields. Throughout this paper (with the exception of the 
enhanced anomalous diffusion discussed in Sec. 2.3) it is as- 


s CEO So 


( pc \* (hj [ZBP 
\GeV / \kpc / \^G/ 


K _L : 


v Ib so So I pc 


Z c B Ib 

2 


7+1 


(16) 


I6 cm -vg()d 0 / pc \3 / Ib y (ZB 
s c \GeV/ \kpc/ \pG 


also noted there that the quasi-linear approximation is not fully 
applicable in this case). 

In order that the propagation of CRs is diffusive the pitch 
angle distribution should be sufficiently isotropic. This is given 
if the scattering length / SC att(p) = v-’/v ;/ is small compared to 
any CR density scale-length (f/(df/dx)). Within the above 


Note that these diffusion coefficients are valid under the as¬ 
sumption that the turbulence dissipation length is smaller than 
the gyro-radius. If this is not the case, a much higher k\\ and a 
much lower k ± would result. 

Cross field diffusion is strongly inhibited. The small mo¬ 
bility of particles perpendicular to the field lines can be am- 
plified by rapid parallel dif f usion along diverging field lines 
( Rechester and Rosenbluth , 1978 ). A particle’s microscopic 
displacement from its original field line of the order of the 
gyro-radius by microscopic cross-field diffusion grows expo¬ 
nentially while it follows its new field line. This can be de¬ 
scribed in terms of a Liapunov length d| , which depends on 
the behaviour of the magnetic field autocorrelation function at 
small displacements, since the old and new field lines will be 
strongly correlated. As soon the particle is sufficiently far from 
its original field line, it can be regarded as de-correlated from 
it. During the time needed for the particle to de-correlate from 
its present field line, the particle is tied to it and follows its 
stochastic wandering while diffusing along it. The stochastic 
field line wandering can be described as a diffusion process, 
with a diffusion coefficient Db which gives the perpendicular 
displacement (squared) per unit length travelled along the field 
line. Thus, during a de-correlation time the particle is doing 
perpendicular steps of the length given by the field line wander¬ 
ing times the typical diffusion length along the field line (dur¬ 
ing this time). Since this leads to a stochastic displacement, 
the combined propagation can be described as a macroscopic 
diffusion process, which is called anomalous diffusions. This 
anomalous diffusion allows the CRs to use the usually more 
powerful large-scale magnetic fluctuations for their transverse 
pr opagation. A math ematical description of this can be found 
in Duffy et al. ( 1995 ), which should be consulted for details. In 
their formalism the anomalous diffusion depends strongly on 
the parameter 


2 As can be found by an order of magnitu de estimations of the pitch 
angle scattering frequencies as given by e.g. Schlickeiser (2002) 


A = 


SB 2 (I b )A || 
V2 B 2 s(r g ) dj. 


(17) 
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where Ty, A ± are the correlation length of the field fluctuations 
along and perpendicular to the main field direction (usually 
/ B ~ /l|| > T ± ). A is therefore mainly the ratio of the mag¬ 
netic turbulence energy density on the largest scales (l B ) to that 
on the small gyro-radius scale. For the assumed magnetic tur¬ 
bulence spectrum (Eq. |l3]) one gets 


A = 


So d|i 


V2e(r„)/lj_ ^/2soA ± \ZeBl B 


pc 


(18) 


The anomalous cross field diffusion coefficient is 


K-d - k± + 


2 D b k\\ 
Tl In A 


(19) 


for time-scales longer than the particle-field de-correlation time 
). D b = 6B 2 (1 b ) A\\/(4 B 2 ) is the field line 
wandering diffusion coefficient, and Tl = B 2 l 2 ± /(A\\ 8B 2 {I b )) is 
the Liapunov length-scale[| This is the amplification length of 
(the expectation value of the square of) a small displacement 
of a particle from its initial field line while travelling along its 
new field line. The characteristic perpendicular length-scale I ± 
of the size of chaotic , non-trivially twisted magnetic field fluc¬ 
tuations is expected by us to be of the order of l B , although 
arguments exist that it might be much smalleijj. With these def¬ 
initions one gets 


(Duffy et al., 1995 


Kd~ k ± 1 + 


A 2 4 

InAZ 2 


( 20 ) 


3 Strictly speaking, only the perpendicular component of the field 
fluctuations SB 2 ± (I b ) enter the equations for D B and A l. However, the 
difference between the perpendicular and total fluctuation levels are 
not very large. Since also the diffusion coefficients depend on the per¬ 
pendicular magnetic fluctuations, one can interpret SB 2 as the perpen¬ 
dicular fluctuations right from the beginning and get a consistent for¬ 
malism. 

4 Pertubatively (to lowest order in SB 2 ) one can show that Z ± is given 
by the second Taylor coefficient of the magnetic field autocorrelation 
function perpendicular to the main field direction: ( 8B(x ) • SB(x + 
£±)) = SB\I b ) (1 - i 4/^1 + 0(4))- More strictly speaking, I ± seems 
to be the Taylor length of the along the mean field direction integrated 
and by dy normalised autocorrelation function (of the chaotic part of 
the fluctuation spectrum). However, if the autocorrelation function is a 
direct product of a parallel and perpendicular profile, then this is iden¬ 
tical with the simplified definition given here. Otherwise it can be ex¬ 
pected to be closetothatlfoweverAthasbeen^ointed^omb^Narayan 
and Medvedev ( [2001 ) that if the chaotic (topologically non-trivially 
twisted) magnetic fluctuations extend from the turbulence energy in¬ 
jection range down to the smallest turbulence scale, the anomalous 
diffusion seems to become very efficient. If all the magnetic power 
is in such chaotic modes the effective cross field diffusion coefficient 
seems to be close to the order of the parallel diffusion coefficient. In 
the case that the magnetic turbulence in radio cocoons is completely of 
this chaotic nature our results would be changed in a way that particle 
escape from radio cocoons is always very fast. Thus the observational 
tests of particle escape mentioned in the introduction may confirm or 
refute such a picture. In the following it is assumed that the small- 
scale turbulence is dominated by travelling Alfvenic waves, which do 
not disturb the topological properties of the main field. Only on the 
largest scales we do expect strong non-trivial topologies, produced by 
the turbulence driving forces. Their strength should be limited (and 
therefore possibly confined from the smaller-scales) by reconnection 
events. Therefore we expect l ± ~ l B . 


v l B So A j! 


pc 


6 eo In A / 2 \Ze BI 
cm 2 v So T 


i-r 


( 21 ) 


: 1.2 • 10 


27 


CEO 


'S. ( p c y i At mv 

z 2 \GeV/ \kpc/ \pG ) ' 


The ratio of the anomalous to parallel diffusion coefficients 




2 In A Z 2 


0.037 S 2 0 Al/l 2 ± 


( 22 ) 


is nearly independent of the particle momentum, and it is in¬ 
dependent of the fudge factor sq. The numerical result is in 
good agreement with detailed Monte-Carlo simulation of par- 
ticl e diffusion in slab-geometry and homogeneous turbulence 
by Siacalone and Jokipii (1999). These authors find a ratio of 
kJk\\ = 0.02...0.04 for a turbulence strength comparable with 
Z>o = 1 (see their Fig. 3) and a scaling which seems to be only 
slightly weaker than oc <5 q (see their Fig. 4). 


2.3. Enhanced anomalous diffusion 

Under certain circumstances the anomalous cross field diffu¬ 
sion is extremely efficient. The anomalous diffusion coefficient 
k li depends quadratically on A (Eq. pf] ), which is mainly the 
ratio of large-scale to small scale turbulence (Eq. |I|). A is 
therefore independent of the level of the turbulence, as long 
as the spectral shape of the turbulence energy spectrum is not 
changed. Increasing the turbulence energy density of a system 
by a factor Xj increases K. d also by X-\ . due to the dependence 
of k_ l oc Xt (or /C|| oc X. y 1 ) on the small scale turbulence. 

If the slope of the turbulence is changed, drastic changes 
in K d can result. If e.g. large-scale kinetic energy is suddenly 
injected into a system, it may need a period on the order of 
the eddy turnover time Tt = Zt/vt before the turbulent cas¬ 
cade has also raised the small-scale turbulence level. During 
this period the microscopic diffusion coefficients k\\ and k ± stay 
unchanged, but /e a is increased by a factor of X 2 . 


3. Cosmic ray escape 


3.1. Boundary layer penetration 

ff a boundary layer of insulating magnetic fields (without any 
flux leaving the boundary layer regions) separates the interior 
from the surrounding of the source region the escape of cosmic 
ray particles is strongly inhibited. Although radio polarization 
observations of radio cocoons indicate that magnetic fields are 


well aligned with the cocoon surface (e.g. Laing, 198C, 1981 
[Spangler et al., 1984), it is far from obvious if this implies an 


insulating boundary layer. The origin of such a hypothetical 
layer may be due to mixing of radio plasma with the ambient 
gas in a thin surface mixing layer, which is then dynamically 
decoupled from the turbulent interior due to its much higher 
inertia. The expected shear between the interior and the mix¬ 
ing layer may amplify and align magnetic fields into such an 
insulating boundary layer. 

As a simplified model we assume the source region to be 
spherical with radius r s , the boundary layer to have a thickness 
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<a?b- The layer, which extends from radius r, - r s - d b to r s , is 
assumed to be filled with magnetic fields which are mainly tan¬ 
gential to the source surface at r s . Therefore the radial diffusion 
coefficient is that of the macroscopic cross field diffusion K d . 
In a quasi-stationary situation (slow particle escape) the space 
density of particles is given by 


g(r, P ) = go(p) 



(23) 


within the boundary layer (r; < r < r s ), where goip) is the 
particle spectrum inside the layer (r < rp. From this the number 
of CRs within the source region is given by 


Nip) 




2 71 


dV s g(r, p ) = — g 0 (p) r s r t (r s + n), 


(24) 



and the CR escape rate through the source surface area A s is 
given by 


Nip) = A s K. d Vg(r,p)\r=r s 


47TK d g Q (p) 


The escape frequency follows to be 


(25) 


Fig. 1. Geometry of the inter-phase region. The left (x < x s ) 
and right regions (xi < x ) are the source and loss regions of the 
particle. A bundle of field lines connect these regions through 
a transition zone (x s < x < x\). 


Vlayer(/?) 


Nip) _ 6 K. d (p) 
Nip) r 2 - r 2 


3 Kaip) 


r s d b ( 1 -d b /(2r s )) 


(26) 


ignored. These conditions lead to the following parameters: 
q(x,p) = H(x s - x)f s ip)/T s ip), (27) 


4.0 v do ^P c 1 / p c y 

(ZB\ 2 

t(x, p) 

= H(Xs - x) Tsip) + H(x - Xl) T\ip) , 

(28) 

Gyr ceq r s d b l 2 ± VGeV/ 

Uoj ' 

p(x, p) 

= 0, 

(29) 


It should be noted that enhanced anomalous diffusion due to 
a sudden raise of the large-scale turbulence level by a factor At 
increases the particle escape rate through a boundary layer by 
a factor X 2 during the period of not fully developed turbulence, 
and At afterwards (see Sect. 2.3). 


3.2. Flux tube escape route 

If some magnetic flux </># leaks from the source to the loss re¬ 
gion, the CR can escape following the field lines. A possible 
situation is sketched in Fig. |T]. The leakage requires that the 
particles of the source region enter the inter-phase flux tube ei¬ 
ther by perpendicular diffusion from disconnected regions, or 
by travelling along the flux tube from a distant reservoir. Then 
they diffuse through a transition zone between the two phases 
and finally leave the tube in the loss region by cross field diffu¬ 
sion (or escape to infinity along the field line). 

For simplicity it is assumed that within the source and the 
loss regions all parameters of Eq. are spatially constant. The 
source function q is given by f s , which is the average density of 
particles within the source (per magnetic flux) divided by the 
time r s a particle needs to cross-field diffuse over the distance 
given by the diameter of the flux tube. The typical time a par¬ 
ticle requires to enter or leave a flux tube with flux </> = n r 1 B 
and radius r ~ /lj_/2 is given by r = r 2 /(4/c a ) = 0/(4 n B K d ). 
In the loss region q — 0. In the transition zone K\\ and B might 
be functions of the position, but no particle losses or sources 
are assumed there. Furthermore, any energy loss processes are 


B(x) = H( x s - x) B s + H(x - x\) B\ 

+H(x - x s ) H(x i - x) B t (x ), and, (30) 

K\\(x,p) = H(x s - x)K % (p) + H(x - X\)K\(p) , 

+Hix - x s ) H(x i - x) *r t (x, p). (31) 


H{x) is the Heaviside step function. Since the momentum is 
conserved (p(x,p) — 0), Eq. |5] can be solved independently 
for all momenta. The explicit dependency on p is dropped in 
the following. The solution of Eq. [s] with the parameters given 
above, smooth field lines ( dB/dx = 0 at x s and jcj) and the 
boundary conditions /(- oo) = / s , /(oo) = 0 is 


(f s ~ ifs ~ /a) exp((x - x s )/l s ) ; x <x s 
f (x)—i [B s f d - ( B s / a - B\ /b) cr(x)] /Bfx) ; x s <x <xi, 
l/b exp((xi - x)/h) ; xj<x 


(32) 


with 


Is — y/Xs Ts , /1 — y/x\ Ti , 

K s I K s l\ B\ 


1 + 1 / 

K\I 


/a = fs/ 
fb = fs/ 

+,) - i 

1 J.r, Xfx'. 


Is Bs + Is B , 


_ | B\ | K\ l s 

I B s B s k s l\ 


and I = 


r;, b^) 

I dx —-— . 

J.v, K t (x') 


(33) 

(34) 

(35) 

(36) 
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Fig. 2. Profiles of f(x), g(x) = B(x) fix), B(x), and Ky(x). The 
parameters are f s = 100, B s = 0.3, B\ = 0.01, x s = -2, x\ = 2, 
Z s = 3, h = 3. The field in the transition zone is parametrised by 
B t (x) = [B\ - (Bi - B s ) e xp(- \[KjB\ Or - x s ) 2 /(x i - x) 2 )] + B s - 
(B s - B\) exp(- \IB\IB S (x - xi) 2 /(x s - x) 2 )/2. Further k\\(x) = 
0.01 B~ 2 ^(x) is assumed. 


Fig. 3. As Fig. j|, but B s = 0.01, B\ = 0.3. 


Strongly magnetised loss region: 

_ %B S pT _ rjs B s k\ _ (p B K\ 

V " ~ L s B\ Vrj " L s l\ B\ ~ BiVJi 


(42) 


An example density profile is shown in Fig. [l] for strong 
fields in the source and weak field in the loss region, and in 
Fig. I for the reversed geometry. 

If one introduces the spatial average within the transition 
zone by (A(x)) = dxA(x)/l t , where Z t = x\ - x s , one can 
write 

I = B t Zt/kt = B t xfrjkt, where (37) 

B t = (B ), K t = B t /(B/K t ), and r, = Z 2 /k t . (38) 


If B changes by a large factor in the transition zone one expects 
max(Zf s , B\) > B, » min(B s , B\). Due to the expected anti¬ 
correlation of B and Ky also max(/i s , k\) » kt > min(/c s , k\) is 
likely. 

The total particle loss of the source into the loss region is 
given by 



_ fs 4 >b B s _ 

Bs Vrs/Ks + Bt xln/kt + B\ VtiM 


(39) 


If several loss channels exist with similar properties (r s , n) then 
4>b in Eq. ^can be replaced by the sum of the absolute fluxes 
of these channels. In this case the flux leaving the source region 
is (/>b = A , Tj s B s , where q s is the fraction of the source surface 
(A s ) penetrated by inter-phase magnetic flux. After defining the 
characteristic length-scale of the source by L s = V s /A s , where 
V s is the source volume, and noticing that N s = V s g s = V s f s B s 
one can write the typical particle escape frequency vy as 


My _ xjs B s IL s 

N s B s y/rJiTs + B t V^t/kt + B\ VnM 


(40) 


This formula is simplified if one of the terms in the denomi¬ 
nator dominates. Two cases correspond to typical astrophysical 
situations: 

Strongly magnetised source region: 


TJs_ _ ?7s k s _ 4>B Ks 

Bs y Ls Is Bs kg Is 


(41) 


If two regions 1 and 2 with volumes Vi and V 2 are con¬ 
nected by some magnetic flux <ps, then the ratio of the particle 
exchange frequencies is equal to the inverse volume ratio, as it 
is required by detailed balance: 

_ ?7i B\!L\ _ V 2 
V2—>1 112B2IL2 Vi 

In the following it is assumed that the source field region 
has the strongest magnetic fields and therefore Eq. [id] is used. 
If a number of My flux-tubes with diameter A ± each leave the 
source region the total magnetic flux leaving the source is <Pb = 
M\\n B T”/4. Applying the identity My = 4 / 7 S A s /(7t A^) one 
finds that the particle loss frequency is then given by 


4/? s *11 

4j_ Ls 

4 Tjs Z B T|| v I pc 

3 V2 In A e 0 /_l dj_ Ls \ZeBl B ) 

0.081 v Z 2/3 Ay /(A ± l ± L s ) I p C /(Z B) | 
~My^ Tls ceo kpc- 4/3 \ GeV//,Ci ) 


(44) 

(45) 


It should be noted that enhanced anomalous diffusion due 
to a sudden raise of the large-scale turbulence level by a factor 
Xj increases the flux tube particle escape rate by a factor X T 
only during the period of not fully developed turbulence, but 
not afterwards (see Sect. 2.3). 

If the CRs escape rapidly from the source region, they can 
excite plasma waves which again scatter the CRs. This limits 
the escape velocity to approximately the Alfven- or thermal ve¬ 


locity, or whichever is larger (Tademaru, 1969; Holman et al. 


1979). For the high energy part of the CR spectrum the lim¬ 


iting velocity can be even larger, since the small number den¬ 
sity of high energy CRs in typical power-law distributions does 
not lead to an efficient excitation of the scattering waves (e.g. 
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Felice and Kulsrucj , 2001). However, an important quantity is 
the CR escape velocity at the surface of the source region: 


i'cr - — = 


4 v 


Ib '*|| 


h 3 so V21n A 4±1± \ZeBls 


pc 


1-7 


(46) 


km v 
80 — 


;5/3 i 

^ D ''Ml 


l pc/(ZB ) y 
S £oc 2 ± /j_ kpc 2/3 \ GeV /pG ) 

As long as this velocity is below the limiting speed, which for 
radio galaxy cocoons should be the environmental sound speed, 
CR self-confinement due to wave excitation should be negligi¬ 
ble. 


3.3. Cross field escape route 

In order to model the cross field escape frequency an individual 
flux tube of diameter A± is investigate first, which touches the 
inter phase surface on a length-scale H|| (see Fig. |4j). The parti¬ 
cle distribution function along this tube is governed by Eq. |5] if 
one adopts the following parameters: 


q(x,p) = (l-cr ± (x))f s (p)/T s (p), 
r(x,p ) = r s (p), 
p(x, p) = 0 , 

B(x) = fi s 
K\\(x,p) = K s (p). 


(47) 

(48) 

(49) 

(50) 


cr ± (x) is the fraction of the flux tubes surface which is in con¬ 
tact with the other phase. cr ± (x ) = 0 for x < xi or x > X2- 
Again, the notation of the momentum dependence is dropped 
for convenience. The solution can be found with the help of the 
method of Green’s functions and is given by 


fix) = fs 


i r Xi 

1 - —— I dx' cr ± (x') exp(-|x - x'\/l s ) 
^ Is JX] 


(51) 


In the following it is assumed that cr ± (x) = const within the 
length A\\ — X2 - x\. In this case 


fix) = fs 


l-2± 

2 


■ xi 


h 


■ h 


X '2 ' 


, with 


lifz) = sgn(z) (1 - e 1:1 ), 


(52) 


which is shown in Fig. |5j. The particle escape rate (into the other 
phase) of this filament is then 


-N ± = 


where 




fs 4|| CTj 


1 -— (y + e y - 1 ) 

>’ 


4 <5 0 


(j. V2 In A 


(53) 


(54) 


Usually y«l, otherwise the situation resembles more the iso¬ 
lating boundary layer case discussed in Sect. 3.1 The number 
M ± of flux tubes which touch that part of the surface of the 
source region, which is not intersected by inter phase magnetic 
flux, is 


M , = 


-4s (1 ~ ds) 

n A ± /Iii cr _l 



Fig. 4. Geometry of a flux tube which is part of the source 
boundary. The left (x < xi) and right regions (x 2 < x) are the 
source and loss regions of the particles. In between, escape into 
the other phase is possible. 


This leads to perpendicular escape frequency of 


Vj. = - 


M ± N ± 

N s 


= Vl,o( 1 - 77s) 


1 - — (y + - 1) 

y 


with 


4 s /lj. _ 4 K- d 
4 V s r s L s A ± 

2 5 0 Ib$v I pc | 

3 £o In A l\ A± L s \ZeBI B ) 

0.016 5 0 v l B 3 Zj/(l 2 ± A ± L S ) l pc /(Z B)f 
~Myr" kpc - 4/3 \GeV/^G/ 


(56) 


(57) 

(58) 


The reverse process, in which particles enter the strongly mag¬ 
netised region, has a frequency which can be calculated from 
the escape frequency (Eq. ^6|) and the condition of detailed bal¬ 
ance (Eq. |43| ). 

It should be noted that enhanced anomalous diffusion due 
to a sudden raise of the large-scale turbulence level by a factor 
Xj increases the cross field particle escape rate by a factor X~ 
during the period of not fully developed turbulence, and Xj 
afterwards (see Sect. 2.3). 


3.4. The best escape route 

The most efficient escape route can be found by comparing the 
escape frequencies for the flux tube and the cross field escape 
route: 

Kl _ /*a~ l -Tfs 1 

V|| V *11 T h 
do 4|| 


cr± 


(y + e- y -l) 


(59) 


(55) 


5/j. 77s 


(60) 
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Fig. 5. Profile of f(x) for cr ± = 0.5, Ty = 2 and different l s as 
labelled in the graph. 


siders further that such radio plasma is buoyant and therefore 


leaves the cluster centre within a few 100 Myr (Churazov et al. 


2001; EnBlin and Heinz, 2002), the fraction of CRs injected 
into a cluster core is expected to be small under such circum¬ 
stances. This can impose constraints on several of the theoreti¬ 
cal considerations listed in Sect. 1.1, which rely on an efficient 
escape of CR electrons, positrons or protons from radio co¬ 
coons into the centre of galaxy clusters. 


Further, during the buoyant rise of a radio cocoon through 
a cluster atmosphere, the cocoon expands adiabatically. If all 
magnetic length-scales scale linearly with the cocoon size L s , 
and the magnetic field strength decreases adiabatically as B oc 
L~ 2 then all escape frequencies discussed here (vi aye r, vy, v ± ) de- 

—2/3 

crease according to v oc L s . This implies that CRs are better 
confined in a larger cocoon, or that CR escape is most rapidly 
at an early stage of the cocoon’s buoyant voyage through the 
intra-cluster medium. 


During the injection phase of turbulence (a few tt) the ratio 
Vj_/vy is increased by an additional factor of X T over the esti¬ 
mates given in Eq. Under certain conditions (bo < 5 tj s < 
Xt bo) the dominant escape mechanism can switch temporarily 
from flux tube escape to cross field escape during the onset of 
turbulence. 

4. CR escape from radio plasma 

4.1. A quiet environment 

The old radio plasma cocoon of a small radio galaxy in the 
central IGM of a galaxy cluster might be of approximately 
spherical shape. Here the following parameters are assumed: 
B = IOjuG, 1b = /ly = 3kpc, l ± = A± = 1 kpc, and a cocoon 
diameter of 50 kpc. The turbulence inside old radio plasma is 
likely the turbulence induced by the environment, and therefore 
assumed to be Kolmogorov-like. The diffusion coefficients of 
10 GeV CR particles are then K\\ = 7 ■ 10 28 cm 2 s _1 /(eo bo), 
k_l = 1.7 • 10 16 cm 2 s _1 eo b 0 , and K a = 1.9 • 10 28 cm 2 s _1 bo/eo- 
If the radio cocoon has a boundary layer of thickness db - 5 
kpc, which is isolating the interior from the exterior due to the 
lack of any magnetic flux exchange, then the escape frequency 
of these CRs is vi aye r = 1-7 Gyr 1 bo/eo- This is extremely slow 
for a low level of magnetic turbulence (bo < 0.01). 

If such a isolating layer does not exist, CR particles can 
reach the surface much easier and leave by cross field diffu¬ 
sion (over a short distance) or by following open field lines. 
Since only reconnection events of the fossil fields with the 
much weaker environmental fields could have opened the mag¬ 
netic topology, 77 s «: 1 is assumed due to the possible rar¬ 
ity of such events. Under such conditions the escape frequen¬ 
cies of 10 GeV particles are vy » 57 GyU 1 rj s /e o, and v± ~ 
31 Gyr 1 bo/eo- The CR streaming velocity along open flux 
tubes is 460 km/s /sq and therefore likely much below the clus¬ 
ter sound speed, so that CR self-confinement is not impor¬ 
tant for these particles. For a low level of magnetic turbulence 
(bo < 0 . 01 ) and a mostly closed field topology ( 77 s < 0 . 01 ) the 
escape of 10 GeV particles requires several Gyr. If one con- 


4.2. Suddenly injected turbulence 

If a cluster merger event suddenly injects large-scale turbulence 
the radio plasma can become transparent even for low energy 
CR particles. The merger might produce turbulent flows with 
velocities of vx ~ 1000 km/s on a scale of /x ~ 100 kpc which 
increase the turbulent magnetic energy density on large-scales 
by a factor of Vx = 10 from initially bo = 0.01 for example. For 
roughly an eddy turnover time Tx = /x/v’x ~ 100 Myr the small 
scale turbulence is not increased in regions far from shock 
waves. During only a tenth of this period (10 Myr, which is of 
the order the turbulent cascades needs to transfer energy from 
the scales of the radio cocoon (10 kpc) down to the CR gyro- 
radius length scales) an enhanced anomalous cross field escape 
should allow roughly 30 %/eo of the 10 GeV CR particles ini¬ 
tially confined in the radio cocoon to escape. A similar num¬ 
ber of the CRs would escape in the remaining 90 Myr of fully 
developed turbulence. Significant losses would occur on even 
shorter timescales in the case that the pitch angle scattering ef¬ 
ficiency is low (eo <sc 1). The total loss of particles escaping 
along inter phase flux tubes is roughly 6 % ( 77 s / 0 . 01 )/e'o <s 1 
during 0.1 Tx, and likely negligible if only a few kpc 2 of the 
cocoons surface is opened ( r <s 0.01). If the cocoon would 
have an isolating boundary layer of thickness 5 kpc, then only 
1.7 % of the CR particles would be released during 0.1 Tx- 


4.3. Merger shock waves 


Cluster mergers happen frequently (Mohr et al., 1995; Jones 
and Forman, 1999| ; Schuecker et al., 2001|) so t hat shock waves 
appear often in clusters (Quilis et al., 1998; [Vliniati et al. 


2000 ), which should raise the turbulence level on all length 


scales during their passage. Therefore the strongly enhanced 
anomalous diffusivity discussed above is not expected to ap¬ 
pear in this case. However, the turbulence level should increase 
substantially, leading to some enhancement of the particle es¬ 
cape rate, and the typical lengthscale L s the particles have to 
travel before leaving the source also decreases substantially. 
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Both effects in combination should lead to efficient particle es¬ 
cape. 

The escape length L s = V s /A s decreases due to two ef¬ 
fects: First, the radio plasma volume shrinks by a factor C = 
V 2 /V 1 = (P\ //F) 3 ' 4 for radio plasma with a relativistic equa¬ 
tion of state which is adiabatically (due to the high internal 
sound speed) compressed in an environmental shock wave with 
pressure jump P 2 IP 1 . Second, the shock wave should disrupt 
the radio plasma into filamentary or torus-like morphologies, as 
seen in numerical simulations (EnBlin and Brugge 11 , 2002) and 
in high resolution radio maps of some cluster radio relics (Slee 
et al., 2001). The numerical simulations show that a spher¬ 
ical radio cocoon transforms into a torus with major radius 
equal to the radius of the original sphere (the torus diameter 
is roughly that of the former sphere). Using the dimensions of 
the torus given by this observation and the above given volume 
shrinking factor C one finds that the characteristic length-scale 
L s = V s /A s decreases due to shock compression and change of 
the morphology by a factor 



(61) 


Thus, for a typical merger shock wave with Mach number of 
the order 3, and pressure jump roughly 10 the characteristic 
length scale decreases by a factor of nearly 6. If at the same 
time the level of turbulence increases by an order of magnitude, 
the cross field diffusion escape rate is enlarged by nearly two 
orders of magnitude. Therefore also for radio plasma passing 
through a shock wave efficient CR particle escape is expected. 


5. Discussion 

A model has been presented describing the escape of CR parti¬ 
cles from regions which are poorly magnetically connected to 
the environment. The main application of this theoryis to co¬ 
coons of radio galaxies. 

In order to apply this model quantitatively to astrophysi- 
cal systems, several poorly known parameters have to be de¬ 
termined, like the geometry of the source ( V s , A s , L s , ;/ s , cr s ), 
the field fluctuation length-scales (d||, A ± , l ± , Ib), the nature and 
the level of the turbulence (<5o), and the coupling strength be¬ 
tween small-scale magnetic turbulence and the gyro-motion of 
charged particles U; (l ). Several of these parameters are closely 
related and can roughly be estimated (e.g. the characteristic tur¬ 
bulence length-scales). Other parameters will remain unknown 
until a truly detailed knowledge of radio plasma is available. 

However, even without detailed knowledge of these param¬ 
eters, insight into the qualitative behaviour of CR transport be¬ 
tween different phases is provided: 

- An isolating magnetic layer at the surface of the CR confin¬ 
ing region can efficiently suppress CR escape compared to 
the case in which flux tubes from the source interior touch 
(or even leave) the source surface. 

- CR escape from a mostly magnetically confined region al¬ 
ways requires parallel and cross field particle transport. 


The escape frequency of cosmic rays travelling along inter¬ 
phase flux tubes is independent of the level of turbulence 
(Eqs. [f4f|45]). 

The cross field escape frequency increases linearly with in¬ 
creasing turbulence level (Eqs. ^6|]58| ). 

The parallel and perpendicular escape frequencies are 
mainly determined by the properties of the strong field re¬ 
gion. 

All escape frequencies considered here have an identical 
scaling with the CR momentum, for any turbulence spectra. 
E.g. for Kolmogorov turbulence it is vy/j_ oc p l/3 , and for 
Kraichnan turbulence it is vy/ ± oc p 1 ^ 2 . 

In the case where the pitch a ngle sc attering frequency is 
strongly reduced ( Chandran . 2000b| ) due to the possible 
anisot ropic nature of weak magneto-hydrodynamic turbu¬ 
lence (Sridhar and Goldreich, 1994: Goldreich and Sridhat, 


1997), this model can still be applicable as long as the parti¬ 


cle transport is diffusive. In such a case one expects so <K 1. 
If the chaotic (topologically non-trivially twisted) part 
of the magnetic fluctuations extends down in the power- 
spectrum to the smallest length-scales, the anomalous cross 
field diffusion may be much larger than estimated here 
(Narayan and Medvedev, 2001) leading to a very rapid CR 
escape from radio cocoons. If this is indeed the case is not 
yet clear, but may be testable by observations of escaped 
particles^} 

In contrast to a statement by Chandran ( 2000 e ) it is shown 
here that the macroscopic diffusive transport of particles 
along a flux bundle with small scale inhomogeneities in the 
field strength is reduced by the presence of strong field bot¬ 
tlenecks (Appendix A). 

A sudden large-scale turbulence injection can lead tem¬ 
porarily to the enhanced anomalous diffusion regime. This 
increases both escape frequencies, but the cross field escape 
frequency by a larger factor. 

Cluster merger events can lead to strong particle losses 
from radio plasma. Shock waves disrupt the relativistic 
plasma into smaller pieces from which CRs can escape 
more easily. It should be noted that cluster merger waves 
and turbulence themself are able to accelerate CRs, which 
can mask (and/or energise) the CRs escaping from radio 
plasma. 

Merger can also lead to efficient particle escape from radio 
plasma not directly affected by shock waves. The enhanced 
anomalous diffusion due to the large-scale turbulence in¬ 
jected by the merger could also lead to efficient particle es¬ 
cape. 


Given the large astrophysical importance of the problem of 
cosmic ray escape from radio plasma cocoons, further studies 
of this problem - theoretically and observationally - would be 
of great value. The presented work provides a starting point for 
such investigations. 
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5 E.g. annihilation of escaped positrons, decaying pions produced 
by relativistic protons, or even simply a low level diffuse synchrotron 
halo around a bright radio cocoon due to escaping electrons. 
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Appendix A: Homogenisation of the diffusion 
equation 


It is instructive to derive the macroscopic transport equation in the 
case that the spatial fluctuations in the parameters of Eq. |I] can be 
described by a large-scale x' and a much smaller scale f, so that 
x = x' + f/e and e is the ratio of the small to the large-scale. 
A quantity A(x ) = A(x J .f) is averaged over the smaller scale by 
<A>(.0 = lim^„ fid?A(x?,?)/£. 

The evolution of the distribution function can be expressed as an 
asymptotic series in the small parameter e. To lowest order in e this 
method of homogenisation gives 


df d - - 


d /_ d(Bf) \ 
dx I* 11 Bdx ) 



(A.1) 


where the prime was dr opped (x' —> V) . The details of this kind of cal¬ 
culation can be found in Holmes ( 1995 1. The homogenised parameters 
are given by 


/ = </> (A. 2) 

B = (l/By' (A3) 

P = B(P/B > (A.4) 

*11 = B/<B/*y> (A.5) 

f = B-'/(1/(tB)) (A.6) 

<? = (?>• (A.7) 


In the case of vanishing small-scale fluctuations of all coefficients, 
their homogenised aver ages are, of course, unchanged by the ho¬ 
mogenisation and Eq. AT is identical to Eq. ||. In the case of 
small scale variations, the homogenised loss coefficients p and f are 
weighted towards the weak field regions, where the flux tube has a 
larger diameter and therefore most of the particles reside. The effective 
large-scale diffusion coefficient ky is reduced compared to the ty pical 
small scale coefficient (/Cy), due to the harmonic weighting in Eq. |A.5j. 
The weight of these regions is further increased due to the expected 
anti-correlation of field strength and diffusion coefficient. However, 
even in the case that K\\ does not change on the micro-scale, its global 
value is changed due to field strength fluctuations on the micro-scale: 
i<\\ = /cy/((l/B>(B» < K\\ (equality only if B(x,f) = B(x)). If e.g. 10 % 
of the micro-scale has a 100 times enhanced field strength, the aver¬ 
aged diffusion coefficient is a factor of 10 smaller than the microscopic 
one, in the case of a spatially constant K\\, otherwise much more. This 
demonstrates that magnetic traps modify the particle transp ort even in 
the diffusive regime, in contrast to an opposite statement by Thandrar| 
(J2000ah. 


Appendix B: Glossary 


A s : source region surface area 

CMB : Cosmic Microwave Background 

CR : cosmic ray 

C : volume shrinking factor of a shock compressed radio cocoon 
So : relative strength of magnetic turbulence at reference scale (Eq. 


D m/j : Fokker-Planck pitch angle diffusion coefficient 
d b : thickness of insulating boundary layer d b = r s - r ; 

So : fudge factor relating turbulence energy density to pitch angle 
scattering efficiency (Eq. |k|) 


s : ratio of scattering frequency to gyro-frequency (Eq. [IJ) 
q s : fraction of A s occupied by inter-phase magnetic flux 
/ : isotropic particle phase-space density per magnetic flux (Eq. 

/ s : average density of particles (/) within the source 
y : spectral index of magnetic turbulence spectrum 
g : isotropic particle phase space density per volume 
K\\ : parallel diffusion coefficient (Eqs. 10 Bn 
k ± : perpendicular diffusion coefficient (Egs/fTf], |l6| ) 

K d : anomalous cross field diffusion coefficient (Eqs. [T^, |2c| ) 

A: Eqs.0.0 

T|l : parallel correlation length of the field fluctuations 
A ± : perpendicular correlation length of the field fluctuations, typ¬ 
ical flux tube diameter 

Ib ■ magnetic turbulence reference scale (Eq. O) 
l± : (8B(x) -6B(x+£]_)) = 8B 2 (l B ) (1 - i^//x + 0(^)) if only the 
chaotic (topologically non-trivial twisted) field fluctuations are re¬ 
garded. 

I s (I\): characteristic distance a CR travels before leaving a flux 
tube in the CR source (loss) region (Eq. |33|) 

L s = V s /A s : characteristic length of source 
4catt : CR pitch-angle scattering length / scatt = v/v^ 

My : number of inter phase flux tubes 
m : particle mass 

: number of flux tubes forming the boundary 
p : cosine of the CR pitch angle 
N s : number of CR in source region 
Vy : CR escape frequency along flux tubes (Eqs. 

Vf, : particle-wave pitch-angle scattering frequency 
: perpendicular escape frequency 
4>b : magnetic flux 
p : particle momentum 

q : isotropic particle phase-space injection rate per magnetic flux 
(Eq.§) 

r g : particle gyro-radius 

r, : inner radius of (spherical) boundary layer 

r s : source radius in spherical approximation 

cr ± : fraction of flux tube surface in contact with loss region 

t : particle escape time from a flux tube 

r s : flux tube escape time in source region 

tj : eddy turn-over time 

v : particle velocity 

14 : source volume 

X T : increment factor of turbulence energy density 
y = Ty/4 : field correlation length over particle diffusion length 
’ (Eq. 0) 

Z : particle charge 
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